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The possible discovery of s± superconducting gaps in the moderately correlated iron-based su- 
perconductors has raised the question of how to properly treat s± gaps in strongly correlated su- 
perconductors. Unlike the d-wave cuprates, the Coulomb repulsion does not vanish by symmetry, 
and a careful treatment is essential. Thus far, only the weak correlation approaches have included 
this Coulomb pseudopotential, so here we introduce a symplectic N treatment of the t — J model 
that incorporates the strong Coulomb repulsion through the complete elimination of on-site pairing. 
Through a proper extension of time-reversal symmetry to the large N limit, symplectic- A/" is the first 
superconducting large N solution of the t — J model. For d-wave superconductors, the previous un- 
controlled mean field solutions are reproduced, while for s± superconductors, the SU{2) constraint 
enforcing single occupancy acts as a pair chemical potential adjusting the location of the gap nodes. 
This adjustment can capture the wide variety of gaps proposed for the iron based superconductors: 
line and point nodes, as well as two different, but related full gaps on different Fermi surfaces. 



I. INTRODUCTION 

The new family of iron-based superconductors^ has ex- 
panded the study of high temperature superconductors 
from the single band, d-wave cuprate superconductors 
to include multi-band superconductors with full gaps. 
Experimental^ and numericaP'^ work suggest a range of 
correlation strengths between different materials. From 
the theoretical point of view, weak and strong corre- 
lation approaches converge on many of the major fea- 
tures: most importantly, the predominantly s± struc- 
ture of the superconducting gapP^^^. The real materials 
are likely in the regime of moderate correlations where 
both approaches are useful. Unlike the d-wave cuprates, 
where the strong Coulomb repulsion is eliminated by 
symmetry, these multi-band s± superconductors require 
a careful treatment of the Coulomb pseudopotential^^. 
While thi s has bee n incorporated into weak coupling 
approache d^ * -'^^ 1 -'^^ ^, it has yet to be included in strong 
correlation treatments based on the t — J model. With 
this in mind, we introduce the symplectic- t — J model. 
The use of a large limit based on the symplectic group, 
SP{N) allows a proper treatment of time-reversal in the 
large- A/" limit^^^^, making this the first superconducting 
large N treatment of the t — J model. Symplectic- also 
replaces the usual U{1) constraint of single-occupancy 
with a SU{2) constraint that strictly eliminates on-site 
pairing. This SU (2) constraint is essential to the treat- 
ment of s± superconductors, where it acts as a pair 
chemical potential, adjusting the gap nodes to eliminate 
the Coulomb repulsion. 

This paper is intended as an introduction to the 
symplectic- A^ t — J model, illustrating the importance of 
the additional constraint with the example of s± super- 
conductors, and showing how this model contains a range 
of gap behaviors reproducing different iron-based materi- 



als. We begin by reviewing the Coulomb pseudopotential 
in section [ll] and demonstrate the lack of superconduc- 
tivity in the SU (N) t — J model corresponding to the 
usual slave boson mean field theory. In section III, we 
introduce symplectic Hubbard operators, which allow us 
to develop a superconducting large- A/" treatment of the 
t — J model. We demonstrate this mean field theory 
on several examples in section |IV| before discussing the 
range of possible future directions in[Vl 



II. THE COULOMB PSEUDOPOTENTIAL AND 
THE t - J MODEL 

On-site pairing is disfavored by the Coulomb pseu- 
dopotential, which will cost a bare energy, UN{0), the 
average of the Coulomb repulsion, V(r^— r^) = e^/|r^— r^l 
over the Fermi sea. However, in the weak coupling limit, 
where we assume the pairing is mediated by the exchange 
of a boson with characteristic frequency ub^ the time 
scale of the pairing, ujb is much longer than that of the 
Coulomb repulsion. In other words, while the effective 
electron-electron interaction is attractive, it is also re- 
tarded, meaning the electrons like to be in the same place, 
but at different times, while the Coulomb repulsion is a 
nearly instantaneous repulsion of two electrons at the 
same place and time. The Coulomb pseudopotential is 
therefore renormalizecP^, 



N{0)U 



l^N{0)U\og 



(1) 



to weak coupling. If Tc cx cj^ exp(— 1/A), the attrac- 
tive interaction is replaced by A ^ A — /i*, which re- 
duces Tc slightly at weak coupling, but does not de- 
stroy superconductivity. In BCS superconductivity, the 
bosons exchanged are phonons, and the Debye frequency, 
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ujd ^ Ep- However, in more strongly correlated super- 
conductors, the two time scales are of the same order, 
and the Coulomb pseudopotential can drastically affect 
the superconductivity. Strongly correlated examples, like 
the cuprate and heavy fermion superconductors, avoid 
this problem by developing a d-wave gap, where the pair- 
ing with a positive gap is exactly cancelled out by that 
with a negative gap, as guaranteed by the d-wave sym- 
metry. This choice of gap neutralizes the Coulomb pseu- 
dopotential. However, the iron-based superconductors 
are widely believed to have an s± gap, where the amount 
of cancellation between positive and negative gap regions 
is not protected by symmetry, and depends strongly on 
the Fermi surfaces. When this cancellation is incomplete, 
/i* reduces Tc and it is extremely important to consider 
this effect when mapping out the phase diagram, as it 
affects the relative stability of s- and d-wave gaps. These 
effects have been incorporated in the weakly correlated 
solutions^ but not yet in the strongly correlated 

approaches. 
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FIG. 1: (Left) In the d-wave gap, the cancellation of the 
superconducting order parameter over the Fermi surface is 
guaranteed by symmetry, as the positive (blue) regions will 
exactly cancel the negative (white) regions. (Right) However, 
in the s± superconducting gap, the amount of cancellation is 
extremely sensitive to the Fermi surface. 

Here, we take the strongly correlated limit, ^ oo to 
eliminate double occupancy, which corresponds to taking 
/i* ^ oo. The Heisenberg model describes the insulating 
half-filled limit of the t — J model, but generally holes 
(n < 1) or electrons (n > 1) will hop around in an anti- 
ferromagnetic background. Doubly occupied states must 
be avoided, and the hopping is not that of free electrons. 
Rath er, it is projected hopping , described by the t — J 



The Hubbard operators, Xah = where \a) = 

|0), ensure that only empty sites, or holes can hop 
(or for n > 1 that electrons can only hop from doubly 



occupied sites to singly occupied sites). Here, X^ro are 
projected hopping operators 

Exact solutions of the t — J model are unavailable, and 
the typical approach is to write down a mean field solu- 
tion using the slave boson approach^^ which divides 
the electron into charged, but spinless holons and neu- 
tral spinous. The most common choice is the U{1) slave 
boson representation^^. 



(3) 



so-called because it is invariant under U{1) gauge trans- 
formations. However, mean field solutions do not neces- 
sarily satisfy all the conditions on the full model, and may 
not maintain the /i* ^ oo limit. Large N approaches 
generate mean-field solutions by extending the SU{2) 
t — J model to some larger group. When developing a 
large treatment of the hopping term, one must take 
care that the two terms are consistent, or in other words, 
that the charge fluctuations described by the t term gen- 
erate the spin fluctuations in the Heisenberg term. The 
algebra of these operators, given by 



(4) 



extends the algebra of Hubbard operators from SU{2) 
to SU{N), and we see that two charge fluctuations in 
sequence give rise to a spin fluctuation described by the 
SU{N) spin operator Xai3 = faf^- 

The large N limit of the fuh SU{N) t - J model is: 

{ij) ij 

Decoupling the J term yields a dispersion for the spinous, 
but not pairing^^. There is no superconductivity in this 
large N limit. 



III. THE SYMPLECTIC-iV t - J MODEL 

A superconducting large N limit requires a proper def- 
inition of time-reversal, as Cooper pairs can only form 
between time-reversed pairs of electrons. The inversion 
of spins under time-reversal is equivalent to symplectic 
symmetry, and the only way to preserve time-reversal in 
the large N limit is to use symplectic spins^^ 



Sap = flfp - al3fLsf- 



(6) 



where a ranges from —N/2 to N/2 and a = sgn(a). 
Here we use the fermionic representation because we are 
interested in the doped spin liquid states that become 
superconductors. Introducing doping means introducing 
a small number of mobile empty states. When an elec- 
tron hops on and off a site, it can flip the spin of the 
site. Mathematically, this implies that the anticommuta- 
tor of two Hubbard operators generates a spin operator. 
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In a symplectic-A/' generalization of the t-J model, anti- 
commuting two such Hubbard operators must generate a 
symplectic spin, satisfying the relations: 



(7) 
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where the last equality follows from the traceless defi- 
nition of the symplectic spin operator, 5*0, ^3 = X^^ — 
-^^a(3' When we represent the Hubbard operators with 
slave bosons, the symplectic projected creation operators 
take the following form^^, 

Xao = flb + af_^a, (8) 

SO that the other two Hubbard operators take the form 



Xal3 — Sap + 

Xoo = b^b 



a^a. 



(9) 



This double slave boson form for Hubbard operators was 
derived by Wen and Lee^^ as a way of extending the local 
SU{2) symmetry of spin to include charge fluctuations. 
In our approach the SU{2) symmetry appears as a con- 
sequence of the time-inversion properties of symplectic 
spins for all even which permits us to carry out a large 
N expansion. The Nambu notation, = (b^^a^) and 
= ifl^^f-a) simplifies the expressions, as Xao — fa^ 
and the hopping term of symplectic-A^ t — J model can 
be written. 



=-e|(^u^i/-+h.c.) (10) 
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Restricting the spin and charge fluctuations to the 
physical subspace requires that we fix the Casimir of the 
Hubbard operators^^. 



C = Sj 



(11) 



where S'^ = \ ^oc(5 ^oci3{j)Sj3a{j)' A detailed calculation 
(see Appendix) shows that 

C = {N/2f - l-^'], 



(12) 



where ^ = (^'f', ^, ^3) is given by 
^3 

Q!>0 



N 

Uf + Hb - Ua - Y 



(13) 



a>0 



In the infinite-^/ limit, the Casimir, C is set to its maxi- 
mal value, and we obtain the constraint =0. Writing 
out the condition that vanishes, we obtain 



a}-aj 



fUj^ = N/2 



alb 



a>0 



(^fj — afjo 



0. 



(14) 



The first equation imposes the constraint on no double 
occupancy. The second terms play the role of a Coulomb 
pair pseudo-potential, forcing the net s-wave wave pair 
amplitude to be zero when superconductivity develops. 
Under the occupancy constraint, there is only a single 
physical empty state, which is 



\0)=(b^ +a^aflji)\n), 



(15) 



for = 2. The physical interpretation of these terms be- 
comes clearer if we pick a particular gauge. Since we only 
have two flavors of bosons and N flavors of fermions, the 
only way the bosons contribute in the large N limit is by 
condensing. As the bosons are condensed at all temper- 
atures, Fermi liquids and superconductors are the only 
possible states; while this situation is clearly unphysical, 
and will be resolved with 1/N corrections, it allows us 
to fix the gauge in a particularly simple way by setting 
a = and condensing only the b bosons, 6^ = Nx/2 
because the bosons carry all the charge in the system. 
The factor of N/2 makes the doping extensive in N. The 
constraint simplifies to. 



JjaJj—a 
^fj — afja 



N{l-x) 
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0, 
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In a mean field theory, these three constraints are en- 
forced by a trio of Lagrange multipliers A = (A+, A_, A3) 
in a constraint term that takes the form 



Hr 



3 



A3 



f jafj' 



N{l-x) 



-A+(a4/;_J + H.c 
(17) 

The first constraint is clearly recognizable as impos- 
ing Luttinger's theorem. This term is present in the 
conventional U{1) slave boson approactP^. The second 
terms impose severe constraints on the pair wavefunc- 
tion when superconductivity develops, implementing the 
infinite Coulomb pseudopotential. For d-wave supercon- 
ductors like the cuprates, which have been the main fo- 
cus of previous t — J model studies, these constraints 
are satisfied automatically, and at the mean-field level, 
there is no difference between the symplectic-A/" limit 
and many of the pre viously considered uncontrolled mean 
field theoriej ^^l^^ l ^^l However, for s± pairing, these addi- 
tional constraints enforce the Coulomb pseudopotential, 
/i* and have a large effect on the stability of s± super- 
conductivity. 

Once the bosons are condensed, and the Heisenberg 
term decoupled, the spinon Hamiltonian is quadratic. 



■T3 + Ui. 



\Xi 



Jij 



(18) 
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where we have introduced the SU (2) matrix notation, 



Ui. 



Xij ^ij 
^ij Xij_ 



(19) 



Xij generates a dispersion for the spinons, while Aij pairs 
them. The fuh Hamiltonian is given hy H -\- He- The 
physical electron, c'^ ^ {h)p + (a)/ will either hop co- 
herently, forming a Fermi liquid when A is zero, or will 
superconduct when A is nonzero. The mean field phase 
diagram is obtained by minimizing the free energy with 
respect to these mean field parameters, Xij ^^^d A^j, 



Xij 

A,, 



^/f\ f. \ 

jy \JiaJJOt/ 



(20) 



and enforcing the constraint on average, fjafja) = 

and (E, ocf]j]_J = 0, where (• • • ) is the ther- 
mal expectation value. 

The Ji — J2 model will have two sets of bond variables, 
Xr] and A^, where r] indicates a link, {ij). We assume 
that xi and X2 are uniform, and allow Ai and A2 to be 
either s-wave or d-wave. When these order parameters 
are Fourier transformed, we find Xk = Xi7ik + X272k = 
2xi(cx + Cy) + 4.X2CxCy and Ak = V is a combi- 

nation of s-wave and d-wave pairing on the nearest and 
next nearest neighbor links. 



extended s 
and we define c^, 



2Als{Cx + Cy) 
2Ald{Cx - Cy) 

2A2s{Cx^y +C 
2A2d{Cx^y 



x-yj 



4A2sCxCy 

: -4A2dSxSy (21) 



COS k^jd^ s^j 



sin krjtt. The full Hamil- 



tonian (including the constraint) has the form. 
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NJVsxXs 
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where is the Fourier transform of t^j, /7k is the Fourier 

^(A++A_). (A2 is unnecessary 



transform of Uij and Ai 



if A is real). This Hamiltonian can be diagonalized, and 
the spinons integrated out to yield the free energy. 



-27VT^log2cosh 



\X,\ 



NJVsxXs 



,(23) 



where cjk V^l^ /^k' «k = A3 - ^ + Xk, and /3k 
A I ■ Ak- Minimizing this free energy leads to the four 
mean field equations. 



dF/dxn 



tanh^ 
2wk 



ak7»7k 



dF/dA^ 
dF/dXs 
dF/dXi 



L 



L 



tanh % ^ , 4 

tanh % 

^ ' ctk - x/2 = 

k ^'^k 

tanh% 
„ ^ /?k 
k 2wk 



: 0. 



(24) 



The first three are identical to those for the U{1) slave 
boson mean field theories^^, but the last enforces the ab- 
sence of s-wave pairing. Ai acts as a pair chemical poten- 
tial adjusting the regions of negative and positive gap. 



IV. SIMPLE EXAMPLES 




J rt 



FIG. 2: (a) The ti - J2 model, (b) The Fermi surface (holes 
shown in red) for the ti — J2 model at intermediate doping. 
In the superconducting state, the gap nodes follow the dashed 
lines, separating regions of positive and negative gap. (c) The 
superconducting transition temperatures for the ti — J2 model 
both with (solid lines) and without the Ai constraint (dashed 
lines) , for s-wave (blue) and d-wave (green) superconductivity, 
d-wave superconductivity is unaffected by the Coulomb repul- 
sion, while the s-wave transition temperature is decreased. 

Now let us see this constraint in action, applied to sev- 
eral simple cases. First, we shall take the simplest lattice 
to exhibit s± pairing: the ti — J2 model shown in Fig. 
[2] (a). Here, only the next-nearest exchange coupling, 
J2 and nearest neighbor hopping, ti are nonzero, which 
leads to a single hole Fermi surface with the potential for 
either d^y or s± pairing. The superconducting transition 
temperatures can be determined by setting A2s/d = 

0+ 

and solving the mean field equations, (24) for Tc. The 
results are shown in Fig. [2|c), where we have calculated 
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the transition temperatures as a function of doping, x 
both with and without the Ai constraint. The d-wave 
transition temperature is unaffected, as dF/dXi = by 
symmetry, but the s-wave Tc is suppressed. Note that 
the two transition temperatures are identical for x = 0. 
Looking at the gap structure. Figure [2|b), we see that 
Ai has adjusted the gap nodes such that there are equal 
amounts of positive and negative gap density of states, 
eliminating the Coulomb repulsion. As there is only one 
Fermi surface in this example, there are necessarily line 
nodes even in the s-wave state. The energetic advantage 
of a fully gapped s-wave Fermi surface is thus lost, so 
that d-wave superconductivity, which requires no costly 
adjustment of the nodes, becomes energetically favorable 
for this lattice. 




X 



FIG. 3: (a) The ti - ts - ts - J2 model, (b) The Fermi 
surface (holes shown in red) for the ti — t2 — — J2 model 
at both low (single Fermi pocket) and intermediate doping 
(two Fermi pockets). The gap nodes of the superconducting 
state are indicated with dashed lines, (c) The doping phase 
diagram for the h —t2 —ts — J2 model, calculated with the Ai 
constraint (solid lines) and without (dashed lines). There is 
a quantum phase transition between d-wave pairing (green) 
and s-wave pairing (blue) as doping increases and the s-wave 
states become fully gapped. 

However, if there are multiple Fermi surfaces, s± su- 
perconductivity can gap out both surfaces with opposite 
signs. If we tune the ti — t2 — ts hoppings, keeping only 
J2, we can obtain such a Fermi surface, with two hole 
pockets, as shown in Figure |3ja,b). Again, we calcu- 
late the s-wave and d-wave transition temperatures in the 
presence of the pseudopotential terms, showing the phase 
diagram in Figure (sjc). Our one-band approach makes 
this difficult, as the size of the pockets shrinks with in- 
creasing doping. The s-wave order parameter has line 



nodes for low doping, which recede to point nodes and 
then vanish as the Fermi surface becomes fully gapped 
at larger dopings, where the s-wave superconductivity is 
more favorable than d-wave, causing a d-wave to s-wave 
quantum phase transition as a function of doping. If we 
had equally balanced hole and electron pockets at zero 
doping, s-wave would likely win out over d-wave at all 
dopings. 



V. DISCUSSION 

This study of the symplectic-A^ t-J model illustrates 
the importance of incorporating the Coulomb pseudopo- 
tential into any strongly correlated treatment of s± su- 
perconductors. The symplectic-A^ scheme provides the 
first mean field solution of the t — J model that is both 
controlled and superconducting. The large-A/" limit is 
identical to previous mean-field studie^^, but contains 
the additional constraint fields X± which enforce the con- 
straint $ = 0. For d-wave pairing, this constraint is inert, 
as the s-wave component of the pairing is zero by symme- 
try, but this constraint plays a very active role for s-wave 
pairing, acting as a pair chemical potential that adjusts 
the gap nodes to eliminate any on-site pairing. As such, 
these models can capture the full variety of gap physics 
proposed in the iron-based superconductors: from line 
nodes to point nodes to two different full gaps that are 
not otherwise expected in a local picture. Properly ac- 
counting for the adjustment of the line nodes is essential 
when comparing the relative energies of d-wave and s- 
wave pairing states. 

However, the large-TV limit suffers from an over- 
abundance of coherence, due to the ubiquity of the bo- 
son condensation. As such, the only phases captured 
here are Fermi liquids and superconductors, and study- 
ing the effects of 1/N corrections is an important future 
direction. This application is especially relevant to the 
cuprates, where there have been many intesting, but un- 
controlled corrections to the mean field theories, reveal- 
ing pseudogap-like phases for med b y pre-formed pairs 
and incoherent metallic region&S^EH. A controlled 1/N 
study of the phase diagram of the t — J model studying 
the differences between s-wave and d-wave pairing should 
be of great interest. 

While the t—J models taken in this paper illustrate the 
basic effect of the Coulomb pseudopotential on strongly 
correlated superconductors, they are but poor approx- 
imations of the real materials, due to the single band 
approximation. A better theory would involve multiple 
orbitals per site coupled by a ferromagnetic Hund's cou- 
pling, — I J//|5'^i • Sj^'i between spins in different orbitals, 
ji ^ 11' on the same site. Current large- A/" techniques can- 
not treat such a ferromagnetic coupling, but future work 
might introduce an uncontrolled mean field parameter or 
take Jh co, which may prove more tractable. 

Interestingly, while the majority of the iron-based su- 
perconductors have at least two electron and hole pock- 



6 



ets, there are a handful of "single band" materials: there [Xoq;,Xoo] = — -^oa- 



are the end members KFe2As2^ and Ki-xFe2-ySe2^^ 
which appear to have only hole^^ or electron pockets^^, 
respectively; and the sing le layer FeSe, which has a 
single electron pocketP^'^. In this local treatment, 
KFe2As2's single hole pocket must lead to a nodal d- 
wave super conduct oi^^, as in the t — J2 example above, 
where the s± transition temperature is always smaller 
than the d-wave temperature. A d-wave gap is strongly 
suggested by recent heat conductivity measurements^. 
On the other hand, Ki_a;Fe2-2ySe2 and single-layer FeSe 
have electron pockets, which can develop node- less d- 
wave order, as originally discussed from the weak cou- 
pling approactP^ESI Including the Coulomb pseudopo- 
tential could again become important in this d-wave sys- 
tem if the tetragonal symmetry were broken. 

Finally, an intriguing open problem in the iron-based 
superconductors is the relationship between the local 
quantum chemistry and the superconducting ordei^^. 
The strong dependence of the superconducting transition 
temperature on the Fe-As angle^^ suggests that there 
might be a more local origin of superconductivity, similar 
to the composite pairs found in heavy fermion materi- 
als described by the two-channel Kondo lattice^^. These 
two origins of s± pairing could then work in tandem to 
raise the superconducting transition temperature^^, and 
as such a future generalization of this work to take into 
account both the local iron chemistry and the staggered 
tetrahedral structure is highly desirable. Such tandem 
pairing might explain the robustness of these supercon- 
ductors to disorder on the magnetic iron sit^^. 
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VI. APPENDIX 



In this section, we show that the operator combination 

C = -Sa/sS/Sa i^Oai^ao] — {^00)'^ (25) 

commutes with the Hubbard operators, where Xqq = 
Xqo + 1. C is therefore the quadratic casimir of the sym- 
plectic supergroup SP(N|1). We also show that 



(27) 



C = (iV/2)2 - 1 - 



(26) 



in the symplectic slave boson representation. 

The Hubbard operators Xoa, Xao and Xqo, together 
with the symplectic spin operators, Sai3 = X^/s — -^Xxx, 
form a closed superalgebra: 

[XoajSf3^] = XojSaf3 -\- O^PX^^Sa^ 
{Xoq;,Xo^} = 



Greek indices indicate spin indices a G [±1/2, ±3/2, ±j] 
where j = N/A and N is even. For simplicity, we use the 
notation a = —a and a = sgn(Qf). The operator 8^(3 = 
Xai3 — j^^ais^xx IS the tracclcss symplectic spin operator, 
while the subsiduary operator, Xqo = Xqq — Xq,q,. 
This graded Lie algebra defines the properties of the gen- 
erators of the symplectic supergroup SP(N|1). This su- 
peralgebra is faithfully reproduced by the slave boson 
representation 



flip 

b^fa 

6^6 + 



■f a^afa 



(28) 



while the spin and subsiduary operator, Xqq are given by 



00 



Xqo ± -^Xaa = ± a'^a ± 1. 



(29) 



By inspection, C contains only rotationally invariant 
combinations of the Hubbard operators and each term 
leaves the number of slave bosons unchanged, so that 
it commutes with Saf3 and Xoo- We now show by di- 
rect evaluation that it also commutes with the fermionic 
Hubbard operators Xq^ and X^o 

First we evaluate the commutator between X^o and 
the spin part of the Casimir, 

[Xoa^ Si3rySji3] = [Xoa, S 13 j]Sji3 -\- S 13 j[Xoa^ Sji3] 

= Xo^S^a ± Xo^S'a/3Sgn(a/3) 

± S(3aXo(3 ^sgn{aj)Sa^Xo^ (30) 

Using the identity Saf3 = — sgn(a/3)5'^^, we can convert 
this expression into the form 

[^Oa^ Sf3-fSjf3] = Xo^Sja ± ^O^^^a + ^/5q!-^0/3 ± S^a^O^ 

= 2{Xo/3, 5'/3a}. (31) 
Next we evaluate 

— [^Oaj [Xo(3,X(3o]] = —Xo(3{XoajX(3o} — {Xoq, , X/3o}Xo/3 
= —{Xo(3,S(3a} — {Xoa.Xoo}. (32) 

Finally, 

[Xoa , -^00 ] = ~ l^Oa , -^Oo]-^00 " -^00 [^Oa , -^00 ] 

= {^Oa^^oo} (33) 



Adding ([31|, ([32| and (|33|) together gives 

[Xoa.C]=0. 



(34) 



Since X^o = Xq^ and C = C"*" is Hermitian, it follows that 
[Xq,o,C] = 0. Thus, C commutes with all Hubbard oper- 
ators, and is thus a Casimir of the supergroup SP(N|1) 
generated by the symplectic operators. 
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To evaluate the Casimir, we insert the slave boson form 
of the Hubbard operators. First, evaluating the spin part, 
we obtain 



where 



-- 2 n/(A/' + 2-n/) -4*J.*/ 
Ea>o"/a/I> while 



(35) 



[^oc^ao] = {N -2nf){nb-na)-N. (36) 

Combining the various terms in the Casimir, we obtain 

C = nf{N ^2-nf) - {N -2nf){nb-na) - N 
- 4(^^ + b'^a){^f + a^6) + 46^aa^6 - (Xoo)^(37) 



By regrouping terms, we obtain 

TV 



C -{Uf ^ rib - Ha 



TV 



TV 



) - 2(n5 - Ha) 



- ^t^ + 46^aa^6- (Xoo)^ 
We now introduce the triad of operators 

^3 rif ^rib-Ua- — 



(38) 



(39) 



where ^] = ^3. Alternatively ^(^1+2^2) = and 
^(^1 — i^2) = ^- The Casimir can then be simplified to 

C = -(^3)' + (n6 + n, + l)2 + (TV/2)2-l 

- (^? + ^^)-(n5 + n, + l)2 (40) 

The terms involving and completely cancel out, 
leaving C= (f )' - 1 - (^^^^ 
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